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1 INTRODUCTION 



ABSTRACT 

Using the halo model, we investigate the cosmological Fisher information in the non- 
linear dark-matter power spectrum about the initial amplitude of linear power. We 
find that there is little information on 'translinear' scales (where the one- and two-halo 
terms are both significant) beyond what is on linear scales, but that additional infor- 
mation is present on small scales, where the one-halo term dominates. This behavior 
agrees with the surprising results that Rimes & Hamilton (2005, 2006) found using 
A'^-body simulations. We argue that the translinear plateau in cumulative information 
arises largely from fluctuations in the numbers of large haloes in a finite volume. This 
implies that more information could be extracted on non-linear scales if the masses of 
the largest haloes in a survey are known. 
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2 METHOD 

The Fisher informa tion about a parameter a Riven a s et of 
data is defined (e.g. iTeemark. Taylor fc Heavenall997^ as 



9^1n£(QJdata) 
9q2 



(1) 



Valuable cosmological information is encoded in the 
large-scale structure of the Universe. Rimes & Hamil- 
ton (2005; 2006, RH) have studied how much cosmolog- 
ical Fisher information is in the non-linear dark-matter 
power spectrum. As a start, they considered the informa- 
tion about the initial amplitude A of linear power. RH 
measured this information from the cov ariance matrix of 
power spectra from A'^-body simulations. Mei ksin fc White! 
^^SIIlJS''^''''"'^^'^'^*-'' Zaldarriaea & Hui ( 199!|, SZH) and 



where C is the likelihood function of a from the data. 

We discuss the information in the (hereafter, implicitly 
non-linear dark-matter) power spectrum P{k) about the log- 
arithm of the initial amplitude A of the linear power spec- 
trum P''"(fc). In this paper, the word 'information' is im- 
plicitly about A, although this framework also allows other 
parameters to be investigated. The information in measure- 
ments P{ki) and P{kj) of the power spectrum is (RH) 



ICoorav fc HtJ (|200ll . CH) have also explored this covari 
ance matrix, without explicitly using the information con- 
cept. RH found that, as expected, information is preserved 
on large, linear scales, but the behavior is surprising on 
smaller scales. On translinear scales {k ~ 0.2 — 0.8 h^^ Mpc), 
the information is mostly degenerate with that in the lin- 
ear regime, but there is significant extra information on 
smaller scales. This does not mean that surveys reaching 
only translinear scales are useless, but that translinear scales 
contribute little extra information about j4 in a survey ex- 
tending into the linear regime. Such clearly demarcated be- 
havior in different regimes is suggestive of the halo model, 
in which the power spectrum consists of a term on linear 
scales and a term from virialized haloes. 
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We approximate the expectation value of the middle term 
on the right side of this equation (i.e. the Fisher matrix) 
with the inverse of the covariance matrix, [C~^]ij. Here, 
Cij = {AP{h)AP{kj)), where AP{ki) is a fluctuation of 
an estimate away from the mean in P{ki). This approxima- 
tion is good if the distribution of estimates of power about 
their mean is Gaussian, which RH showed to be sufficiently 
so in this case. 

We follow the procedure of RH to measure I a- We use 
uncorrelated band powers B{ki) with a diagonal Fisher ma- 
trix. The band powers are decorrelated in such a way that 
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they retain the same expectation values as P{ki). Following 
RH, we use an upper-Cholesky decomposition, ensuring that 
the information up to a given fcmax is affected only by data 
for k ^ fcmax- The cumulative information Ia{^ k) which 
we measure in terms of band powers is then 



Ia(^ k) 



Y^ 91n5(fc) d^\uC dlnBjk) 
2-^ dlnA d\nB{kY d\nA 



(3) 



2.1 Covariance matrix construction 

The covariance of the power spectrum in a survey of vol- 
ume V is the sum of a Gaussian term, which depends on 
the square of the power spectrum itself, and a term involv- 
ing the (hereafter, implicitly non-linea r) trispectrum (SZH; 
iHamilton. Rimes fc ScoccimarrolboOfti : 
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where Va,i is the volume of shell i in Fourier space (propor- 
tional to kf for logarithmically spaced bins), and Tij is the 
trispectrum averaged over shells i and j; 

cPk cPk ■ 

T[ki,—ki,kj,—kj)— — . (5) 



K,. Vs. 



To obtain the power spectrum and trispectrum, we use 
the halo-model trispectrum formalism developed by CH, 
which we sketch here. In the halo model, the matter power 
spectrum is the sum of one- and two-halo terms. 



P(fc) 



:P'''(fc)+P^''(fc) 



=M2°(fe,fe)+p""(fc)[Mll(fc)]^ 

where Mj^ are integrals over the halo distribution; 



(6) 

(7) 



M^iki 
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h(3{m)n{m,c) 
xu{ki,m, c) ■ ■ ■ u(ki_i, m, c) dcdm. 



(8) 



Here, p is the mean matter density, m is the halo 
mass, c is the halo concentrat ion as in the NFW profile 
jNavarro. Frenk fc Whitelll99d), baim) is the /3-order halo 
bias JMo. Jing fc Whitelll997l:IScoccimarro et al.ll2001^ . and 
u{k,m,c) is the halo profile in Fourier space, normalized 
to unity at fc = 0. Like CH, we use a distribution in halo 
concentration since its effect can be large for higher-order 
correlations. 

The halo-model trispectrum is the sum of one-, two-, 
three-, and four-halo terms (CR: ICooravll200ll . appendix). 
We quote those terms which seem to have a simple physical 
explanation (Sect. l^TTl : the one-halo term. 



1 \rvi^ Ki^ Kj ^ Ki J — ik/4 t^rCi , fCj , fCj , h^j J 



(9) 



(the arguments of T terms will be implicit henceforth), and 
the component of the two-halo term which comes from tak- 
ing three points in one halo and one point in the second. 



T31 = P'''\h)Ml{h, kj,kj)Ml(h) + 3 perm. 



(10) 



In implementing the halo-model trispectrum, we use the 
same cosmological parameters as RH: (Om, f^A, /b, ft, fs) = 
(0.29,0.71,0.16,0.71,0.97). The only halo-model parameter 
change from CH in our fiducial model is our use of the 
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Figure 1. The Fisher information in the halo-model dark- matter 
power spectrum up to wavcnumber fc, as defined in Eq. |3] for dif- 
ferent variants of the halo model. In each panel, the dotted curves 
show what RH found from hundreds of Af-body simulations of box 
sizes 128 and 256 h~^ Mpc. (For consistency, for the 128 h~^ Mpc 
box size, we use only the results of their 200 particle-mesh sim- 
ulations.) Like RH, we normalize all /a(^ fc) for a 256 h~^ Mpc 
box. The solid black lines are from the fiducial model, with all 
halo-model inputs the sam e as in CH excep t for th e mass func- 
tion, which we take from ISheth fc TormenI ll999); CH used a 
|Press-Schechter 1 1974) mass function. The top panel shows that 
the mass function controls the shape of the translinear plateau. 
The dashed green curve is from an ST mass function calculated 

using P''°(fc < ,„„ , "^-1 ,, — ) = 0, truncated to mimic the con- 

^ ^ 128 h ^ Mpc ' 

ditions in a finite box. Changes in the small-scale slope of halo 
density profiles (middle panel) or in the slope 7 of the relationship 
between halo mass and concentration parameter (bottom panel) 
alter the small-scale shape of Ia{^ fc)- lu the bottom panel, the 
curves rise as 7 grows more negative. 



ISheth fc TormenI ljl999l ST), instead of the lPress-Schechteil 
l|1974) . mass function. We evaluate the trispectrum in the 
center of each bin, performing only an angular average. For 
the fiducial model, we tried decreasing the bin size by a fac- 
tor of five, which barely affected the cumulative information. 
Using all the same parameters as CH, we obtained ex- 
cellent agreement with the power spectrum they plotted. 
We also verified that our code gives the correct relevant 
M^, functions in an a nalytical version of the halo model 
iCoorav fc Shethll2002l . §4.3). We used the convenient form 
of the angular-averaged perturbation-theory trispectrum 
provided by SZH (their eq. 7), and checked our perturbation- 
theory trispectrum against simple results they gave. Also, 
our r'^'' agrees with the covariance in P^^ due to halo 
Poisson sampling, as theoretically expected (see appendix). 
However, our square-configuration trispectrum Aaq(fc) was 
slightly lower than that plotted by CH; also, our correla- 



Halo-model dark-matter power spectrum information content 3 



tion matrix '^ appears to reach a given value at a 

slightly larger k than in Table 1 of CH, by a factor of ~ 1.5. 
But these discrepancies are noticeable in the regime where 
the simple, well-tested T term dominates. After extensive 
tests, we are confident that our code is accurate. In any case, 
our qualitative result is insensitive to these slight differences. 



3 RESULTS 

Figure Q shows the cumulative information 7a (^ k) as 
found with the halo model under different sets of param- 
eters, and as measured by RH from A'^-body simulations. 
All curves share these qualitative features: on linear scales, 
7a (^ k) oc k'^, as expected; on translinear scales, there is 
little additional information; and on the smallest scales, the 
information rises again, but less steeply than in the linear 
regime. This small-scale rise seems to be shallower in the 
halo model than in simulations, although the small-scale 
measurements by RH were near their resolution limit. 

Generally, the halo mass function affects the details 
of the translinear plateau in the cumulative information 
7a (^ k), and halo density profiles affect the small-scale up- 
turn. Changing the inner slo pe of the halo den sity profile 
between -1 (NFW) and -1.5 jMoore et all200ll) . fixing the 
density at the scale radius Va, gives little change in 7a (^ k). 
We changed the slope 7 in the relation c(m) = co{m/mt)~' 
between halo mass and concentration parameter more dras- 
tically; vai2ingj]^_^|om_zero to four times its fiducial value, 
—0.13 JBullock et al.l200H) . has a large effect on 7a(^ A:) on 
small scales. When c{m) does not change with m, there is 
an apparent plateau in 7a (^ k) at the smallest scales, while 
7a (^ k) is greatest when c{m) falls off most sharply. There- 
fore, it could be said that the information on small scales is 
encoded in the systematic variation in halo density profiles 
with halo mass. 

An interesting question which RH posed is whether the 
information (e.g., about A) in the power spectrum is pre- 
served in time. If it is, then the cumulative information mea- 
sured up to a fixed large physical (not comoving) wavenum- 
ber should remain constant in time. It would not make sense 
if information were created, but non-linear evolution could 
decrease with time; for example, it could divert into higher- 
order statistics. Requiring information to be preserved or 
lost is a test of the model which could be used to constrain 
halo-model parameters. 

Our fiducial model passes this test; Figure |5| shows in- 
formation curves as a function of physical k for the fiducial 
model at fairly low redshifts (where the halo model is well- 
tested), along with curves from RH for comparison. Look- 
ing at the right edge of the figure, it seems that the halo 
model predicts that information up to a fixed physical scale 
of 100 hMpc~^ gradually decreases with time. 



3.1 Physical explanation: fluctuations in the 
number of large haloes in a finite volume 

One reason to investigate the information concept with an 
analytic model is to gain physical insight. Much of the be- 
havior of the cumulative information seems to come from 
the following Poisson halo abundance fluctuation (PHAF) 
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Figure 2. The information content of the halo-model power 
spectrum about the initial amplitude of the linear power spec- 
trum, as a function of physical scale, for various redshifts. Also 
shown are the results from 200 particle-mesh simulations by RH, 
for boxes of size 128 h~^ Mpc. For a = 1, we also show the re- 
sult if the linear power spectrum used for the ST mass function 
is truncated to mimic a box of this size. The information up to 
the physical scale of 100 h Mpc"'^ in this fiducial model decreases 
slowly with time. 



effect: in a finite volume V, the number of haloes of mass 
between m and m -I- dm is not Vn{m)dm (which would give 
fractional haloes), but is an integer drawn from a distribu- 
tion with mean Vn{rri)dm. We assume that the distribution 
is adequately Poissonian for haloes large enough that only 
a few may occur in the volume, and that the large number 
of smaller haloes is stable enough that non-Poissonianity in 
the distribution has a negligible effect. 

In fact, the one-halo term T^^ of the halo-model covari- 
ance corresponds to the covariance ( AP'^'^(fci)AP^'^(fcj)\ in 
P^^ from Poisson-sampling the number of haloes in each 
bin. This can be shown analytically by considering mass 
bins small enough to contain only zero or one halo (see ap- 
pendix). We also verified this numerically by comparing T^'^ 
to the covariance of P^^ for 262144 Poisson realizations of 
the mass function. Compared to the true CH covariance in 
eq. l2J, we consider the PHAF covariance to be an approxi- 
mation which is not entirely accurate, but which shares two 
terms (see appendix) with the true covariance. 

Figure |3] shows the power spectra from twenty Pois- 
son realizations of the ST mass function in a volume V — 
(128 h~^ Mpc)^. (Fluctuating the mass function also causes 
P^^ to fluctuate on large scales, but we have renormalized 
P^^ to match P'"' in the largest-scale bin.) This pleisiosaur 
shape helps explain the shape of the information curve in- 
tuitively. Fluctuations in the number of the largest, rarest 
haloes in a finite volume lead to large variances (and co- 
variances) in P^^ on large scales. On translinear scales, this 
effect squelches the information. On larger, linear scales, in- 
formation is preserved since P^^ dominates, washing out the 
fluctuations in Pih. On fully non-linear scales, the power 
spectrum is dominated by power from small haloes, which 
are more stable in their numbers, giving less covariance, and 
therefore significant information. 

This pleisiosaur shape can also be understood as a dis- 
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Figure 3. Power spectra from twenty realizations of an ST halo 
mass function for a volume V = (128 h~^ Mpc)'^. The number of 
haloes in each mass bin &(m) is drawn from a Poisson distribu- 
tion of mean Vn{m)S{m). The fluctuated 1- and 2-halo terms, and 
their sum, are shown in red, yellow, and green, respectively; the 
unfluctuated power spectrum appears as white dots. (The 2-halo 
terms are renormalized to agree with the linear power spectrum 
in the largest-scale bin.) There is a great deal of variance, and 
covariance, in the translinear regime, since it is here that Pois- 
son fluctuations in the largest (rarest) haloes have the greatest 
influence on the power spectrum. 



persion in where P(k) turns away fr om P''"(fc). This sur- 
gests that the HKLM 'scaling Ansatz' PH amilton et alJl991i 
IPeacock fc Dodds: 1996) may still be valid conceptually, de- 
spite RH's apparent evidence to the contrary, if the dis- 
persion in the function taking linear to non-linear power is 
properly taken into account. 

Figure 2] shows how large a role T^^ plays in the halo- 
model covariance. Keeping the Gaussian part from both 
terms of P{k) fixed, but using only T^^ for T, the infor- 
mation I^{^ k) shares some features with the full Ia(^ k). 
On small scales, T^'^ dominates T, and it is still significant 
on translinear scales. But this is not the whole story, since 
the translinear plateau of Ia^{^ k) is indistinct and wig- 
gly. We have not determined the source of the wiggles, but 
odd behavior in I^{4i k) is not necessarily worrisome, since 
it uses only part of the physical T, yet uses the Gaussian 
covariance from the full P{k). 

There is also a term (P^^(ki)P'^^(kj)\ in the halo abun- 
dance fluctuation covariance (see appendix) which equals 
Til (eq. IIUII . However, it is less certain that T^i may be 
attributed to PHAF. We measured the PHAF covariance 
of the full power spectrum, and found off-diagonal corre- 
lations near unity in the linear regime, which do not oc- 
cur in linear theory. These high correlations come from 
(P^ {ki)P^^{kj)j, in which mass function fluctuations cause 
P^^ to be biased relative to P''° by a constant (in k) on 
large scales. (This bias was removed in Fig. |3) The error 
likely comes from replacing the bias b2{m) in the (true) CH 
covariance with [6i(fc)]^ in the PHAF covariance. This re- 
placement does not occur in T^i, so we still attribute T^i to 
PHAF, but less confidently than we do T^^. Figure 0] shows 
what happens if T|^ is included or excluded along with T^^; 
these two terms dominate the behavior of /a(^ k). 



Figure 4. The information content in the halo-model power 
spectrum about the initial amplitude of linear power, if only T^'^, 
and Tgf along with T^^, are included in the trispectrum. Also 
shown is the information content if all halo-model trispectrum 
terms except these are included. The Gaussian piece of the co- 
variance, using the full halo-model power spectrum, is always in- 
cluded. When covariance matrix terms are added, Ia{^ k) de- 
creases, but not simply additively, since the Fisher matrix is a 
matrix inverse of a sum, not a sum of matrix inverses. The Ih 
term can be attributed to Poisson fluctuations in the halo mass 
function, as can, perhaps, T^^. These two terms dominate the 
behavior of the information curve, indicating that Poisson halo 
fluctuations may largely be responsible for the loss of information 
in the non-linear regime. 



4 DISCUSSION 

The qualitative agreement we found between power spec- 
trum information content (about the initial amplitude A of 
linear power) in the halo model and in A^'-body simulations 
(RH) is a success for the halo model paradigm describing 
the power spectrum and trispectrum. Our results also lend 
weight to the argument made by RH that the small-scale 
upturn in cumulative information, seen in their simulations, 
has a real, physical origin and is not the result of insuf- 
ficient resolution. Howe ver, the agreement is not per fect, 
as with the bispectrum (iFosalba. Pan fc Szapudill2005ri . so 
more halo-model in gredients, such as halo triaxiality (e.g. 
ISmith. Watts fc Sheth 20 06). may be necessary for the halo 
model to reach high precision for higher-order statistics. On 
the other hand, the high (co)variance on translinear scales 
suggests that inaccuracies in the halo model on these scales 
have less statistical significance than they seem. 

Even if our information curves do not match exactly 
what RH found, the halo model allows us to understand their 
results better. The halo model will also allow the information 
about other cosmological parameters than A to be studied 
easily. With a covariance matrix such as we used here, all 
that is needed to investigate another parameter a is — ^j^, 
which the halo model can provide self-consistently. 

Fluctuations in the number of haloes of a given mass 
(particularly, massive and rare ones) in a finite volume ap- 
pear to be largely responsible for the paucity of cumulative 
information on translinear scales. Therefore, prospects for 
extracting cosmological information from the power spec- 
trum on these scales might improve substantially with prior 
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knowledge of the mass spectrum in a survey. For example, 
a conditional power spectrum depending on the mass of the 
largest cluster in a survey could contain significantly more 
information on non-linear scales than the power spectrum 
alone. 
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APPENDIX A: COVARIANCE FROM HALO 
MASS FUNCTION FLUCTUATIONS 

Often in the halo model, integrals over halo mass functions 
are multiplied together. Formally, Poisson-sampling the halo 
mass function introduces an extra term in the expectation 
value of these products. Consider the ensemble average 



f{m)h{m) dm / g{m)h{m) dm 



(Al) 



over fluctuations n{m) of the mean halo number density 
n{m), where f{m) and g{m) are functions (e.g. halo density 
profiles) of mass. A Poisson realization of the mass spec- 
trum may be considered, in a way reminiscent of Peebles' 
(1980) derivation of power spectrum shot noise, by parti- 
tioning the mass function into bins 5{mi) small enough that 
in a given volume V, the number of haloes in each bin, 
(/}i = Vh{mi)5{mi), is or 1. Expression lAH becomes 



V2 



y^ fi9j4>i4>i ) = Y2 (^ fi93<t>^ 



I »7^J 



+ V^W^f^9^'t'^)^ 



since (jyl ~ 4'i- Returning to integral notation, this is 



/ n dm I gn dm + r— j fgn dm. 



(A2) 



(A3) 



Products of integrals over the mass function abound in 
the covariance (APiAPj) of the halo- model power spectrum 
(eq. |7|l, where AP{k) — P{k) — P{k). Considering the prod- 
uct of two fluctuations in P^^, 



(AP''^(fcOAP''^(fcj)) = ^Ml{h,h,kj,kj) 



V 



(A4) 



which happens to be the one-halo term of the full CH co- 
variance, eq. @. Extending the notation of eq. © to allow 
more than one factor his{m) in the integrand (e.g. M2^(fc, fc) 
would have two b\{k) factors in the integrand), the lh-2h 
term is, to order 1/V (there is also a 1/V^ term), 

(AP^'^(fcOAP''^(fc,)) = ^P''-^{k,)[Ml\ki,h)M^{k,,kj) + 

2Ml{ki)Ml{ki,kj,kj)]. (A5) 

The second term here (when combined with the same term 
of {AP^'^{kj)AP'^^{ki))) is the CH covariance from T^^; see 
eq. (Ell. To order 1/V, the 2h-2h term is 



{APf^AP. 



2h\ 



p''''(ki)p'''^(kA 



[M^\ki,ki)Ml{k2f + 



2Mi'{ki,k2)Ml{ki)Ml{k2) + 

Ml\k2,k2)Ml{k2f]. (A6) 

If M2" is replaced with Mi, the M^^ terms in eqs. HA5fl and 
<A6l l correspond roughly to T22 and the part of T^*^ which 
does not involve the bispectrum (see CH). 

Taken at face value, eqs. IjAll \A?i\ have signiflcant im- 
plications not only for the power spectrum covariance, but 
for the power spectrum itself, since the two-halo term in- 
volves the square of an integral over the mass function (eq. 
|7|l. By averaging together P^^ from 262144 Poisson-sampled 
mass functions, we verified numerically that 



{Ml(kf) = Ml{kf + ^Ml\k,k) 



V 



(A7) 



However, the 1/V halo discreteness term cannot con- 
tribute to the two- halo term. The two-halo term is the convo- 
lution of the halo-halo power spectrum with density profiles 
of two separate haloes, but in the halo discreteness calcu- 
lation, there can only be (zero or) one halo in each bin, so 
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the 'J = j' term does not contribute. A similar argument ap- 
plies to terms in the halo-model trispectrum; the only terms 
which contain products of integrals over the mass function 
are multiple-halo terms. The covariance from halo discrete- 
ness, on the other hand, is meaningful because it involves a 
product of complete power spectrum terms. 



